Survival probability and first-passage-time statistics of a Wiener process 
driven by an exponential time-dependent drift 
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The survival probability and the first-passage-time statistics are important quantities in dilTerent 
fields. The Wiener process is the simplest stochastic process with continuous variables, and im- 
portant results can be explicitly found from it. The presence of a constant drift does not modify 
its simplicity; however, when the process has a time-dependent component the analysis becomes 
difficult. In this work we analyze the statistical properties of the Wiener process with an absorbing 
boundary, under the effect of an exponential time-dependent drift. Based on the backward Fokker- 
Planck formalism we set the time-inhomogeneous equation and conditions that rule the diffusion of 
the corresponding survival probability. We propose as the solution an expansion series in terms of 
the intensity of the exponential drift, resulting in a set of recurrence equations. We explicitly solve 
the expansion up to second order and comment on higher-order solutions. The first-passage-time 
density function arises naturally from the survival probability and preserves the proposed expansion. 
Explicit results, related properties, and limit behaviors are analyzed and extensively compared to 
numerical simulations. 



I. INTRODUCTION 

Since the primeval discovery of the Brownian motion 
and its statistical analysis, the variety of applications in 
which a relevant stochastic description result is obtained 
is steadily increasing P, [2| • The Wiener process and the 
Ornstein-Uhlenbeck process are idealized statistical de- 
scriptions that apply to many systems. One of the most 
valuable theoretical tools to deal with these and other 
stochastic processes is the Fokker-Planck (FP) formal- 
ism [1, 0] • In this framework different realizations of a 
system are described by the probability density to find 
the system in a given state at a certain time, and a dif- 
fusion equation describes its temporal evolution. Two 
related questions of wide interest in several areas are 
the probability that the system remains in a certain do- 
main at a given time and the instant at which the system 
leaves it for the first time. Given the stochastic nature 
of the process, different realizations of the system leave 
this survival domain at different times and it is natural 
to consider what the statistical properties of this random 
variable are. This question constitutes the so-called first- 
passage-time (EFT) problem [JH-lll- 

The survival probability as well as the EFT problem is 
easy to formulate but difficult to tackle, except for some 
simple cases. In particular, for a Wiener process driven 
by a constant and positive drift toward a fixed positive 
boundan^, these quantities have a simple analytical solu- 
tion However, the extension to a time-dependent 
drift is not straightforward, mainly because the system 
is no longer time-homogeneous. In one-dimensional sys- 
tems, the main work on this topic possibly is [9]. In 
that study, the author extended, via the forward EF de- 
scription, the classical results of Siegert [1] for a particle 
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being driven by a small time-dependent potential, super- 
imposed on a general field. By applying a perturbation 
scheme, the author derives the recurrence relations be- 
tween the linear corrections of the moments of the EFT 
density function. Other series of works have analyzed 
the behavior of the system in a time-dependent sinu- 
soidal drift, in general, studied in the context of stochas- 
tic resonance (see ^10- 13] for seminal works on this topic 
for the Wiener process with an absorbing boundary; for 
other processes, we refer the reader to 1J|). However, we 
are interested in the EFT problem of the Wiener process 
driven by an exponential time-dependent drift because 
it naturally arises in neuroscience, when modeling spik- 
ing neurons with adaptation currents [isj . This process 
can also be used to model a neuron with an exponential 
time-dependent threshold (see [l6l | for the transformation 
between an exponential time-dependent drift to an ex- 
ponential time-dependent threshold). With reference to 
moving thresholds, the main related work is p7| . where 
the authors analyzed the moments of the EFT density 
function of a Markov process with a moving barrier, giv- 
ing some specific examples applied to biological sciences. 

In this work, we study the survival probability and 
the EFT density function of the described system in the 
framework of the backward EF formalism. We describe 
the complete statistics, instead of focusing on its mo- 
ments as in previous studies. We obtain the equation 
and conditions governing the survival probability and 
propose a solution in terms of an expansion in the ex- 
ponential drift intensity. This results in an infinite set 
of recurrence equations, which we explicitly solve up to 
second order. Higher-order terms are outlined and dis- 
cussed. In particular, we show that all order functions 
exist and depend exclusively on the actual time difference 
when the initial conditions are imposed for the backward 
state, as physical considerations require. This constitutes 
the exact solution of the problem. From the knowledge 
of the survival probability it is straightforward to derive 
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the complete FPT density function, which in turn results 
in an expansion series. Since it is natural to solve the 
equations via a Laplace transformation, we review some 
important related properties easy to compute from the 
Laplace transform of the FPT density function. 

In the second part of this work, we focus on the explicit 
results we have obtained and compare them with numer- 
ical simulations. Given that truncation of the expansion 
series results in an approximate solution, we mostly ana- 
lyze the system in the linear regime. Related properties 
and the behavior of the linear solution in different limits 
are also considered. 



II. THEORY 

In one dimension, the FPT problem can be basically 
formulated as follows: a state variable evolves stochasti- 
cally according to a given law in its phase space, and we 
are interested in describing when this variable leaves a 
certain domain for the first time. To deal with this prob- 
lem a number of different methods or approaches had 
been described, mostly based on the knowledge of the 
time-dependent probability density or its temporal evo- 
lution [3-01 ■ Here, we first solve the survival probability 
in terms of the backward FP equation and then derive 
the FPT density. 



F{t'\xo,to) 



P{x\t'\xo,to) dx'. 



(2) 



To describe the transition probability density we use 
the backward FP equation. In this case, given that the 
state variable is at position x' at time t', the probability 
density of the particle being at the position x at an earlier 
time t (t < t') is given by 



dP{x',t'\x,t) 



dt 



fi+ — e-(*-*")/''<i 



dP{x',t'\x,t) 
dx 



-D 



d^P{x',t'\x,t) 
dx"^ 
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The drift coefficient quantifies the first moment of the 
differential transition density in the neighborhood of the 
backward state {x^t). Necessarily, the local level of the 
exponential term is relative to the initial time tg, break- 
ing up time homogeneity. 

The FPT is incorporated with the initial and bound- 
ary conditions: P{x' ,t'\x,t = <') = 1 for x < Xthr and 
for X > a;thr, and P{x\t'\x = Xthr,t) — 0. 

Integration of Eq. Q in x' from — oo to xthr yields the 
survival probability from time t to time t': 



A. Survival probability 

The nonautonomous system we address here can be 
described in terms of the Langevin equation. 



dt Td 



(1) 



where x is the state variable (position, voltage, etc.; here- 
after, the position), t is the time, fi is the constant part 
of the drift, e quantifies the strength of the exponential 
time-dependent drift with time constant Td, and ^(t) is 
a Gaussian white noise characterized by {^{t)) = and 
(^(0C(^')> = 2D S{t - t'), with Z? as a constant. 

Suppose we have a particle at position xq at time 
and it evolves to a position x' at a posterior time t' 
[t' > to) according to a transition probability density 
P{x' ,t'\xo,tQ). Glearly, in a FPT problem a certain re- 
gion of the domain is forbidden and actually the transi- 
tion probability density has implicitly incorporated this 
fact. In this work, we analyze the region defined by a 
constant boundary, x' < Xthr, which set the survival do- 
main, so the forbidden region is a; > Xthr- 

The survival probability F(t'\xo,to) is the probability 
that the particle remains in the survival domain at time 
t' given the initial conditions, and it is given simply by 
integration of the transition probability in the x' domain 



dF{t'\x,t) 
dt 



/i H e 

Td 



-(t-t0)/Td 



dF{t'\x,t) 
dx 
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d^F{t'\x,t) 
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Since t' is a parameter, we make the substitution 
T — t' — t and rename the probability F{x,T]t'). The 
corresponding equation is 



dF{x,T;t') 
d^ 



dF{x,T;t') 
dx 



-D 



d'^F{x,T;t') 
dx"^ 



(5) 



with F{x, T = 0; t') — 1 for x < xthr and for x > Xthr, 
and F{x = a;thr,T;t') — 0. 

To solve this equation we propose an expansion in pow- 
ers of e: 



F{x,T;t') = Fo{x,T;t') + eFi{x,T;t') + e^F2ix,T;t') + . . . 

OO 

= ^e"K(x,T;t'). (6) 

n=0 

Replacing Eq. ^ into Eq. ([SJ and grouping in orders 
of e, we obtain 
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n=l 



OK 
dr 



dFo dFo ^d^Fp 

dFn 



dx 



dx 



-D 



dx 



.y,2 



= 0, 



(7) 



where we have simpUficd the notation for the sake of 
clarity. 

Since e is a parameter, each term in brackets should be 
identically 0. Therefore, to find the survival probability 
we have to solve 



dFo 
dr 


dFo 
^ dx 


^d^Fo ^ 
dx"^ 


dF^ 
dr 


dFn 
^ dx 


d^F 
D " = 
dx'^ 
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— e 


'(t'-t0)/-!-d 




Td 





= 



dF„ 



dx 



for n > 1. 



(8) 



(9) 



This system of Eqs. ^ and can be solved recur- 
sively up to the degree of accuracy needed. To complete 
the solution of the survival probability we have to define 
the initial and boundary conditions for all the functions 
Fn{x,T]t'), for n > 0. Again, given the arbitrariness of 
e, the nonhomogeneous conditions should be imposed to 
the zeroth-order function. Therefore, initial conditions 
are 



Fo(a:,T = 0;0 



1 if a; < Xthr 
if a; > Zthr 
Fn{x,T ^ 0;t') = for n > 1, 



(10) 
(11) 



whereas boundary condition is Fn{x — xthr, t] t') = for 
n > 0. 

Next we solve the expansion up to the second-order 
term and analyze higher orders. 



1. Zeroth order solution 



where we denote Foi-'^) Laplace transform of Fo{x, r) 
to the s domain (due to time homogeneity, t' only appears 
in r). In deriving Eq. we have used the fact that 
Fo{x — )■ — oo) is bounded. 

By the inverse Laplace transformation of Eq. (fT2|) . we 
obtain the solution in terms of r = — Fo{x, t' — t). At 
this point we state the initial conditions of the problem, 
X — xo and t — to- Therefore, F{x, t' — t) F{xo,t' — to). 
Again, replacing t ~ t' — to (now, r is the actual time 
difference) and transforming back to the s domain, we 
obtain 



/ N 1 1 r 

'(s) = cxp<^ 

s s L 



■ (a:t hr - Xo) 
2D 



}■ 

(13) 

where now we have recovered the dependence on s in the 
notation. By comparing Eqs. and we note just a 
single change, x xo- However, the procedure described 
is important in time-inhomogeneous problems and it will 
be important when solving the following orders. 

Even when the inverse Laplace transform of Eq. (I13p is 
available, we disregard this step since, as we will sec later 
when deriving the FPT density function, it is unneces- 
sary (and actually it is related to the FPT cumulative 
distribution) . 



2. First order solution 



The first order term is given by the solution of Eq. 

for n = 1, 



dFi{x,T;t') _ ^ dFi{x,T;t') _ ^ d^F^ (x, r; t') 



1 

Td 



dx 

-{t'-to)/Td 



d_ 

dx 



dx^ 

e"/^^ Foix,T;t')\ , (14) 



with the corresponding initial and boundary conditions. 

This equation can be solved via the Laplace transform. 
In the s domain, Eq. ([Ti)) reads 



The system described by Eqs. ([8]) and ((T0| corresponds 
to the constant drift case (e = 0). The survival probabil- 
ity of the Wiener process with constant drift and diffusion 
coefficients is a time-homogeneous process (the system 
remains unchanged with a shift in t') and easy to solve 
in Laplace domain. Omitting the dependence in s (to 
solve the equation, s acts as a parameter) , this probabil- 
ity reads 



Fo'^ix) 



■ exp ■ 



(a;tiir 



2D 



H - ^/Jj^ + ADs 



(12) 



s F^{x;t') - 



dFtix;t') _^d^Fl-{x;t') 



dx dx^ 



T"d 



where represents the Laplace transform opera- 

tor and Fi{x; t') is the Laplace transform of Fi{x, r; t'). 
Substituting the result we obtained before, Eq. (|12p (note 
that the initial state of the problem is not already eval- 
uated), into Eq. we have 
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dx dx^ 



i_ Q-{t'~to)/ri 
Td 2D{S - 1/Td 

(Xthr - 2:) 



• exp • 



2D 



V/i^ + AD{s - 1/Td)] }.(16) 



The general solution to this equation is given by 



F^{x]t') = Ci exp - 



H+ y/fi^ + 4Ds 



2D 



/ M - yZ/i^ + iDs 

L2 exp :rT^ X 



2D 



1 



-(t'-t0)/Td 



2D (s ~ 1/Td) 

■ {Xthr - X) 



■ exp ■ 



2D 



fi - + ^D{s - 1/Td)] }, (17) 



vahd for Re(s) > l/Td- 

Taking into account that Fi{x — > — oo;t') is bounded 
and the boundary condition is F^{x = Xthr\t') — 0, we 
obtain 



Ft{x-t') = ,^o-(''-'°^'^^ 



G^-i/Td) 



exp 



exp ■ 



• {xthr - X) 

2D 

• (xthr - X) 

2D 



^Ai^ + AD{s - 1/Td)] } 



+ 4i:'s] } 



(18) 



Fi(a;,T;t') = ^ e-H'~io)/r^ ^ r/r^ 



2ttj 



^y? + \Dz - n 



(7z — J 00 

{xthr - 2:) 



exp 



exp ■ 



2D 

• (xthr - x) 

2D 



/i - + 4^?(Z + 1/Td)] } 

} [ (20) 



Vjj^ + 4D2 



where now, it is easy to check that the region of conver- 
gence of the integrand is Re(z) = tJz > 0. However, there 
still is an exponential factor that makes the expression 
diverge. 

At this point we are able to evaluate the real condi- 
tions of the problem: x = xq and t — t^. Remembering 
that T = t' — t, the latter condition imposes that the 
two exponential factors (before the integral) in Eq. ([20)) 
cancel each other. Hereafter, we use t to represent the 
actual time referred to the initial time, t = t' — t^. Tak- 
ing the Laplace transform on this variable, from Eq. ([20)1 
the function -F'i(t) (note that x was evaluated and the 
dependence on t' is exclusively given by the combination 
in t) transforms to 



1 



2D s 

(xthr - 2^0) 



•{ 2D 

2D 

valid for Re(s) > 0. 



- + ^D{s + 1/Td)] } 



We further operate with the inverse Laplace transform 
of Eq. (UHl), which is 



Fi(x,T;0-^e-(*'-*'')/-- 



- + 4£'(s - 1/Td) 



27rj7<T-joo (s-l/T-d) 
(a;tiir - x) 



exp ■ 



exp ■ 



2D 

{xthr - X) 

2D 



Vm^ + AD{s - 1/Td)] } 
y^^^TlD^] }| ds, (19) 



where j represents the imaginary unit and a > 1/Td. 
Taking the substitution z = s — 1/Td, we obtain 



3. A note on the higher-order solutions 

In this subsection we remark on some aspects of the 
existence and the convergence of higher-order terms ex- 
pressed in the Laplace domain. The higher-order terms 
in the expansion, Eq. ([6]), correspond to the solution of 
Eq. ^ with the appropriate initial and boundary con- 
ditions (n > 2). In particular, we obtain an equation 
analogous to Eq. (|T4| : 



dFn{x,T;t') ^ dFn{x,T;t') ^d^Fn{x,T-t') 



dr dx 

i_ e-(*'-*o)/^<i — 
Td dx 



dx^ 

fl^^ F„_i(x,T;i')] -(22) 



The term on the right-hand side of the equation cor- 
responds to a source because Fn-i{x,T;t') was already 
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solved. As in the first-order case, the knowledge of the 
source term in the s domain enables us to readily Laplace 
transform the equation, obtaining an ordinary differen- 
tial equation with a forcing term. The homogeneous part 
of the solution is exactly the same as that in Eq. p7)) 
(terms with unknown constants Ci) and the particular 
solution is different for different orders. Moreover, Ci 
has to be for bounded solutions and the existence of the 
particular solution is given as a sum of the exponential 
factors present in F^_-f^{x;t')\g_i/T.^. After evaluation of 
the boundary condition, the solution is given as the sum 
of n -|- 1 exponential terms. 

Here we note the structure that this forcing term 
imposes on the solution. Since Eq. ([22|) operates in 
the backward state {x,t), the previous-order solution 
Fn-i{x,T;t') must not be evaluated in the initial state 
{xo,to), or correspondingly, its Laplace transform should 
be done for the variable t = t' — t and not for t = t' — to. 

To simplify, we exemplify the concepts of convergence 
with the second-order solution and then extend the con- 
clusion to all orders. In the Laplace domain, the equation 
governing the second-order solution is 



~r, dF^(x;t') (fFi-{x;t') 

sF^{x;t')-n — ^ : ' ' - D- ' 



dx 



:;^e-(*'-«)/-^{/:[Fi(.,r;0](.)}j._,/.,,(23) 

where the Laplace transform of the previous-order so- 
lution is given by Eq. ((T5| . Prior to the evaluation of 
the initial state, this solution has a region of convergence 
Re(s) > 1/rd. It is easy to check that, due to the de- 
lay introduced in the Laplace domain, the forcing term 
in Eq. ([25)1 will impose that the region of convergence 
of the Laplace transform of the second-order solution 
is Re(s) > 2/rd, but now a factor exp[— 2(t' — to)/Td] 
appears. Explicitly, the equation governing the second- 
order is given by 



~r, ,^ dF^(x:t') d^F^{x:t') 
s F^{x;t') - fi- 2 V 1 ^ n 2 V , ; 



dx 



- D- 



dx'^ 




2D{s - 2/Td) 



• exp 



• exp ■ 



(Xthr - X) 

2D 



- Vm' + - 2/Td)] } i ,(24) 



i^2^^;i')-2^e-^(*'-*°)/- 



M - vV + 4D(7^27^ 



2D{s ~ 2/Td) 



Oi exp I 



i=0 

where 

flo = 

fli = 

02 = 



(a;thr - x) 
2D 



v/m' + 4Z?(s - V^d)] }, (25) 



+ ^D{s - l/Td) - 2 V^m' + - 2/^d) 



Ai- v//-t'+4Z?(s-l/Td), 

1 r 



M - + AD{s - 2/Td) 



(26) 



Proceeding as in Eqs. (ITO)) and (PU)) , we obtain two 
exponential factors, exp[— 2(i' — io)/Td] and exp(2T/Td), 
that cancel each other when the initial state (xo,to) is im- 
posed. This cancellation actually means that the second- 
order term of the survival probability, with the initial 
state imposed, depends on time exclusively through the 
combination t = t' — Iq (actual time difference). There- 
fore, its Laplace transform (on the variable t ^ t' — tg) 
with the initial state evaluated is 



-y^JJ+lDs 



•exp I 
where 

bii-s) 



{Xthi- - Xq) 

2D 



- Vm^T4:D(J+77^] }, (27) 



-V? 



4Ds 



^ti^ + AD{s + l/Td), 
^/li^TiDi^Tl/^ - ^V^i^+4Ds 

-i/i. (28) 



and its solution is 



Recursively, the nth-order solution in the backward 
state {x,t) has a Laplace transform valid for Re(s) > 
n/Td with a factor exp[— n(t' — io)/''"d]- Therefore, the 
preceding conclusion extends to all orders. 

Given the existence of the solution of all terms, the 
expansion proposed in Eq. ([6]) constitutes the exact so- 
lution of the system. 



B. First- passage time density 

In the previous subsection, we demonstrate that the 
expansion given by Eq. ©, with the initial state evalu- 
ated, {xo,to), actually reads 
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F{xo,t) = Fo{xQ,T) + eFi{xo,T) + e^F2{xQ,T) + ... 

OO 

= ^e"F„(xo,r), (29) 

n=0 

where the dependence on time appears exclusively 
through the combination t = t' — tQ. 

Once the initial conditions are stated, by definition, 
F(xo.,t) is the probability that the particle remains at 
time T = t' — Iq in the survival domain and, hence, 
equals the probability that the FPT is posterior to r: 
F{xo,t) = Prob(r > r), where T represents the FPT. 
In terms of the cumulative distribution function of the 
FPT random variable, $(r), this means that F{xq,t) — 
1 — $(t) (hereafter, xq is a parameter and can be dis- 
regarded from notation). The density function, 0(r), is 
given by 



d$(r) 
dr 



dF{xo,T) 
dr 



(30) 



which means that the FPT density function has an ex- 
pansion given by 



0n(T) 



dFn{xo,T) 

dr 



(35) 



In the s domain, Eq. (|34|) reads 



OO 

n=0 



(36) 



where <j>n(s) is the Laplace transform of the nth term in 
the expansion of (/'(t), c/jnir), and it is given by 

4>^{s) = 1-s/'o'^(^o,s), (37) 
4>f,{s) = -sF^{xo,s), forn>l. (38) 

For example, from the findings in the previous subsec- 
tion, the terms in the expansion up to the first order of 
the FPT density function are 



h (s) = exp| — 



^l- ^/ii^TAIX^Y (39) 



n=0 



(31) 



Remembering that the initial condition in the diffusion 
problem reads F{x, r = 0; i') = 1 for a; < Xthi, it results 
that F{xo,T = 0) = 1 in the solution already evaluated 
with the conditions of the problem (and obviously xo < 
Xthr for a nontrivial problem). Therefore, the Laplace 
transform of Eq. ((30|) reads 



(32) 



0^(s) = l-sF^(xo,s), 



where <j)^{s) [F^{xq^ s)] is the Laplace transform of (/)(t) 
[nxo.r)]. 

Equivalently, in terms of the expansion for F^ [xq , s) 
[see Eq. ([29]) ] , the Laplace transform of the density is 



n=0 



(33) 



where Fj^{xo, s) is the Laplace transform of the nth term 
in the expansion of F{xo, r), F„(xo, r). 

Since 0(r) has an expansion given by Eq. pip , it is 
natural to write 



(r) = £e" Mr), 



(34) 



where 



(.) = 



2D 



r (a^thr - Xq) 

2D 



exp 



M - + 4i?(s + 1/7-d)] } 
{^^^^[M-v/^^^TiK]}}. (40) 



Eq. (1391) is the classical result for the FPT problem 
with constant drift /i and diffusion D coefficients 7], con- 
sistent with our approach. 



C. Related properties of the first-passage time 
density 

Since the solution of the proposed expansion is nat- 
urally obtained in the Laplace domain, here we review 
some properties easy to calculate from this knowledge. 
It is easy to check that the moments of the density func- 
tion satisfy 



which means that all the moments preserve the expansion 
in e 



(42) 



ri,=0 
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where 



(43) 



For example, the first two raoments for the unper- 
turbed case {n = 0) are 



2D{xt\,^ - xo) (a;thr - a;o)^ 



(44) 



which correspond to the constant drift case 

The linear changes in these properties, Eq. for 
n = 1, are 



exp 



2D 



- 1 



M(2;thr 



exp 



+ 4i?/Td 
(Xthr - Xa) 



(Sthr 




(45) 



These results, Eq. coincide with those of the 

corresponding case in 9J, obtained from a different ap- 
proach. 

The assessment of the complete density function in the 
Laplace domain enables us to obtain another important 
property. The successive ordering of FPTs of Wiener pro- 
cesses, each of them independent of the history (in our 
system, this means the fixed escape domain and initial 
state), constitutes a renewal process. Given the times 
{tk\ when the system reaches the threshold xtiu- starting 
from xq and setting it again to xqj we can construct a 
"spike train", X{t), defined by 



(46) 



representing a renewal point process. 

The Fourier transform of X{t) is the spike train spec- 
tral density Siio), which represents an important prop- 
erty in some fields, such as neuroscience 22|. It is re- 
lated to the density function of a single escape process, 
expressed in the Laplace domain, through [2l|, [23| 



S{u:) 



2tt(t) 



0^(j^) 



(47) 



where uj is the angular frequency. The sequence of the 
renewal times is mean subtracted; otherwise, a S peak 
appears at frequency 0. 

In particular, the constant driving case (e ~ 0) is an- 
alytically tractable and a relatively simple expression 
is found in [22|. The exponential driving case (e ^ 
0) corresponds to the spike train produced by a per- 
fect integrate-and-fire neuron with an exponential time- 
dependent threshold 



III. COMPARISON TO NUMERICAL RESULTS 

In this section we test different theoretical results and 
compare them with numerical simulations. As shown, 
the expansion given by Eq. ([6]) is the exact solution of 
the system. However, the explicit computation of suc- 
cessive terms in the expansion is performed up to certain 
order. Since truncation introduces an error for any fi- 
nite order, we mainly focus on the first-order expansion 
with small values of e. In this case, the time-dependent 
exponential drift can be thought of as a perturbation to 
the unperturbed system defined by e = (constant drift 
case). Without mathematical loss, we set all quantities 
of the system to nondimensional units. 



A. Linear order expansion 

In Figs. [Ija) and [Ijb) we show the FPT density ob- 
tained from simulations for different intensities of the ex- 
ponential drift [(a) e = -0.5 and (b) e = -2.0]. The his- 
togram obtained from simulations (stair-like solid line) 
is compared with different predictions. The zeroth-order 
prediction, (j)o{T), is given by the inverse Laplace trans- 
form of Eq. (PI 



Mr) 



(Xthr - X q) 



■ exp ■ 



[(Xthr - Xq) ~ Ht] 

4Dt 



■}, (48) 



and corresponds to the constant drift case (e = 0). The 
linear-order solution is composed of (poir) + c'/'i (''')• The 
function (jfi (r) is obtained from numerical Laplace inver- 
sion of Eq. (PU)) . 

As shown in Fig. [Ha), for low intensities of the expo- 
nential drift, the FPT statistics is well characterized by 
the linear order. As expected, when the intensity is in- 
creased, higher-order effects become significant and the 
linear expansion is not enough [Fig. [Ijb)]. In this case, 
the second-order solution, given by Eqs. (P71) and ([55)) 
and numerical inverse Laplace transformation, success- 
fully accounts for the numerical data. 
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From numerical simulations 

— Zeroth order expression 

— First order expression 
Second order expression 



(a) 




2 3 4 

First-passage time x 

FIG. 1: (Color online) First-passage time density func- 
tions </>(t) for different intensities of the exponential time- 
dependent drift. Histograms obtained from simulations (stair- 
like solid lines) are compared with analytical results: zeroth- 
order [dotted-dashed (blue) lines], (j)Q{T); first-order [dashed 
(red) lines], 0o(t") + e0i(r); and second-order [dotted (cyan) 
lines], (jioiT) + e<^i(T") + 6^(^2(1"), expressions. 4>o{t) is given 
by Eq. (|48|l. whereas </'i(r) and <^2(t) are obtained from nu- 
merical Laplace inversion of the corresponding expressions 
[Eq. (gOl) for 01 (r) and Eqs. ^ and ^ for <?!>2(r)]. (a) Low 
intensity, e = —0.5; (b) high intensity, e = —2.0. Insets: The 
analytical expression [solid (green) line] for the first-(second-) 
order function </>i(r) [02 (t)] is compared with the empirical 
linear (second-order) function (stair- like solid lines), see text 
for definition. Remaining parameters are N = lO'^ simula- 
tions for each case, jj, — 1.0, D = 0.1, a^thr ~ xq — 1.0, and 
Td = 10.0. 



To directly compare the linear correction 0i (r) with 
its numerical equivalent, we construct an empirical lin- 
ear function as follows. From N FPT processes for the 
constant (e = 0) and time-dependent (e 7^ 0) cases, we 
obtained their histograms, 0const(''') and ^timcdcp. (t) , and 
define (t) = [(/)timedop.(T)-0const(T)]/e. Obviously, 

^const (''') coincides with (/>o(t) (not shown). If the FPT 
process is dominated by the linear regime, the empiri- 
cal function so obtained should agree with the analyt- 
ical result, (I>i{t). As shown in the inset in Fig. [IJa), 
both functions coincide for a small perturbation e — —0.5 
[stair-like line represents the empirical function, whereas 
the solid (green) line is the analytical expression] and are 
a relative mismatch for e = —2.0 (not shown). In this 
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FIG. 2: (Color online) The first two moments as a function 
of the (negative) intensity of the exponential perturbation, 
— e, for different Td. Analytical expressions [solid (red) line] 
are given by Eqs. H44p and H45p . whereas data obtained from 
simulations are represented by different symbols. Remaining 
parameters as in Fig. [T] 



case, given the empirical linear correction constructed 
with a tiny perturbation e = —0.1, we can construct 
a similar empirical second-order correction function as 

(/)2"'P"'(r) = [(/'timodcp.(T)-(/'const(T)-e (/i^^"' (t )] /e^ . As 

shown in the inset in Fig.[ljb) this empirical function co- 
incides with its analytic counterpart ^lij) (fluctuations 
due to a finite number of simulations become higher than 
in the linear construction). 

Given the parameters fi = 1 and (a;thr — xq) = 1.0, 
the mean FPT for the unperturbed case is 1. As shown 
in the inset in Fig. [Ija), a small additive exponential 
drift generates a biphasic correction to the constant case 
density, with a positive weight for times shorter than 
the unperturbed mean and a negative weight for longer 
times. The overall effect is to distort the density, increas- 
ing the probability of times shorter than the typical time 
in the unperturbed system leaving the survival domain, 
for positive perturbations. For negative perturbations, 
the contrary is true. 



B. Moments of the first-passage time density 

As expressed by Eq. (|42p . a linear correction to the 
FPT density of the unperturbed system is reproduced in 
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000 



|°gio 

FIG. 3: (Color online) The first two moments as a function 
of the exponential time scale rd (semi-logarithmic plot), for 
different strengths of the perturbation. Remaining data as in 
Fig.H 



all its moments. In Fig. [2] we show the first two moments 
for the unperturbed case (e — 0) and different strengths 
of the perturbation, e 7^ 0. The analytical expressions 
for these moments are given by (r'^) = (t'^)^^ + tij^)^^ 
(k = 1,2), where (t'')^, are given by Eqs. ([4l|) and (|i5| . 
Arbitrarily, we use negative values for e. In this case, 
moments shift toward larger values in comparison to the 
unperturbed case. As shown, both moments [Figs. [H^a) 
and m^b)] coincide with the linear analytical results for 
low intensities and mismatch for larger values. To de- 
scribe these properties properly for large values of e, we 
should include higher-order terms in the expansion given, 
Eq. (H^ . The range of validity of the linear regime is 
given by the time scale of the exponential perturbation 
Td, which set the effective intensity in the Langevin equa- 
tion [see Eq. ([Ij]. As the time scale of the exponential 
drift increases, the linear regime remains valid over a 
wider range in the e coordinate. 

To stress the preceding paragraph, in Fig. [3] we show 
the first two moments as a function of the time scale (log- 
arithmic scale) for different perturbation intensities. In 
this case, a given perturbation e produces a linear dis- 
tortion for a large time scale Td, but higher-order effects 
become important for smaller time constants. However, 
as we see in the last subsection, the results we obtained 
for the linear expansion still hold in the limit of vanishing 
time scales Td — t- 0. 



C. Spectral density 

To simplify the following discussion, here we set time 
units in milliseconds, while maintaining nondi- 
mensional magnitude. In this case, D, and Td are 
measured as 1 /milliseconds, 1/milliseconds, and millisec- 
onds, respectively. As indicated in Sec. IlICi the spectral 
density of a renewal process composed of consecutive first 
passages (hereafter called the spike train) is easily com- 
puted with the Laplace transform of the density function, 
Eq. (H71) . For exaniple, the (one-sided) spectral density 
(per unitary time) [23| of the unperturbed system has a 
relatively simple expression (see Eq. (3.17) in [l^]), which 
is shown in Fig.SJ^a). 

As indicated at the beginning of this section, we re- 
strict ourselves to consideration of the effect of a small 
additive exponential time-dependent perturbation on the 
spectral properties of the spike train evoked by a sys- 
tem driven by a leading constant drift fi. For such a 
situation, the change in the spectral density of the un- 
perturbed system is hardly noticeable, and therefore, to 
analyze the frequency-dependent changes introduced by 
the perturbation, we focus on the ratio between spec- 
tral densities. If the unperturbed system is character- 
ized by /i, D and Xthr — [FigHJ^a)], and an expo- 
nential drift is added to the system, defined by e and 
Td, the effect of this perturbation on the spectral den- 
sity is shown in Fig. Sfb). In this case, we set a neg- 
ative perturbation e < 0, which implies that the mean 
FPT increases [see Eq. P5|)]. Consequently, the rate of 
the consecutive first-passage processes decreases to lower 
frequencies, in comparison with the unperturbed case. 
This rate is strictly given by the value of the spectral 
density at infinite (two distant events are uncorrelated, 
which means a "white" spectrum at frequencies tending 
to infinity, for stationary processes), S{lo ^ cjo), and 
it roughly determines the position of the observed peak 
[see Fig. 111(a)]. For example, in Fig. SJ^a) the rate is 
defined by (t)~^ — fi/{xthr — xq) = 1 ms~^. This re- 
sults in an asymptotic value of 10"^ s^^ [a factor equal to 
2 ap pea rs when considering the one-sided spectral den- 
sity [23], such as that shown in Fig. HJa)] and a peak 
located near 10'^ Hz. The lower rate obtained by the ad- 
dition of a negative exponential perturbation decreases 
the asymptotic spectral value and shifts the peak to a 
lower frequency. In the spectral ratio we consider in 
Fig. Ub), /Sport (w)/S'unpcrt(w), thcsc cffccts are reflected 
by an asymptotic value less than unity and a biphasic 
shape, with a positive (negative) peak located at a lower 
(higher) rate than the unperturbed rate. 

The effects already mentioned (representative of the 
linear regime) are mainly related to the mean (t) , so we 
consider an alternative situation where this property does 
not change. In this case, if the unperturbed system has a 
given mean [in Fig.|4l[a), (t) = 1 ms], the perturbed sys- 
tem will be driven by the low exponential time-dependent 
drift, freely defined by e and Td, and the constant compo- 
nent will be modified (in comparison to the unperturbed 



10 



(a) 



a 

o. 

CO 
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FIG. 4: (Color online) (One-sided) Spectral density (per uni- 
tary time) of spike trains for different cases. In each case, 
the spectral density was obtained as an average of 10® inde- 
pendent simulations, (a) Spectral density of the unperturbed 
case (e = 0). Parameters: fi — 1.0 ms~^, D = 0.1 ms~^, 
Xthr — xo ~ 1.0. (b) Ratio between spectral densities, 
Spert((^)/5'unport(tij), for the exponential time-dependent drift 
as a perturbation. Same parameters as in (a) for each case, 
and the perturbed system is additionally defined by e = —0.1 
and Td = 10.0 ms. (c) Ratio between spectral densities, 
<5'pert,mc(w)/Sunpert('»j), for the meau corrected exponential 
drift (see text). Same parameters as in (a) for the unper- 
turbed case. For the exponential drift case, D and a^thr — xq 
are the same as in the unperturbed case, and the exponential 
drift is characterized by e = —0.1 and Td = 10.0 ms [same 
as in (b)], but the constant component of the drift, /ipcrt.mc 
is changed in order to obtain the same mean (r) as in the 
unperturbed case ((r) = 1.0 ms). 



system) in order to keep the mean unchanged. In con- 
sequence, the asymptotic value for the ratio between the 



spectra should be equal to unity, as shown in Fig. |4l^c) . 
In this figure we note that other effects are present since 
the ratio is not flat, but they are an order of magnitude 
less than the case where the exponential term is a direct 
additive effect to the unperturbed system [Fig. |4l(b)]. 
Moreover, these effects are hardly noticeable by simula- 
tions (even for the large set we used). Therefore, the 
change in the mean is the main effect introduced by this 
kind of time-dependent perturbation. 



D. Limit behaviors 

Finally, here we analyze the behavior of the first-order 
analytical solution we explicitly obtained, for large and 
small time scales of the exponential drift. 



1. Limit Trf — > oo 

In section |lT] we obtained the FPT density function as 
an expansion in e and explicitly derived the first-order 
expression [for the FTP density, it is easy to obtain the 
second-order expression from Eq. (|27|) via Eq. (|38l) ]. In 
particular, we have this solution fully characterized in 
the Laplace domain. To analyze the limit of large time 
scales, Td — >■ oo, we expand Eq. (|^ in terms of (1/Td) 
and keep the lowest terms. Up to order 1, this expansion 
reads 



■ exp ■ 



2DTd 



ADs 



2D 



}.(49) 



Prom this expression we can obtain, by differentiation 
in s and evaluation at s = 0, the linear term in the ex- 
pansion of the moments of the density function, valid for 
this limit. The resulting expressions coincide with those 
in i- 

Eq. (jini) shows the Laplace transform of the linear term 
appearing in the e expansion, for the limit Td — oo. In 
particular, it is possible to find its inverse Laplace trans- 
form [201, which reads 



(a^thr - Xq) (xthr - Xq) - /XT 



■ exp • 



2Dt^ 



[{xthr - Xo) - AiT]' 



ADt 



(xthr - Xq) - /iT 

2Dta 



Mr)- 



(50) 



In Fig. [SJa) we show the product of the asymptotic 
linear correction function, Eq. (|50p , multiplied by Td as a 
function of r. This product does not depend on Td and, 
therefore, can be compared on the same scale with the 
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products obtained from simulations. As expected, as 
increases, both the empirical product and the analytical 
result coincide. 

From Eq. ([50]) . it is easy to see that the FPT density 
function in this limit, up to order 1, is 



1 



(xthr - a;o) - HT 



2DTd 



Mr). (51) 



In fact, instead of working out the limit rd ^ oo in 
Eq. ((40|l as we did before, we can look at the actual 
physical situation. For Td — ^ oo, the exponential drift 
can be thought of as a constant (i.e., Td » (t)). In this 
case, the constant drift would be /i + e/rd and the FPT 
density function would be given by (j)Q (r) , Eq. (1481) , with 
this modified drift: 



(r) = /— exp . 



4Dt 



(52) 

Expanding Eq. ([52]) around e = up to order 1 gives 
the same result as obtained before, Eq. ((5T|) . 



(a) 




T^= 1.0 
T^= 10.0 

Analytical expression 



First-passage time i 




First-passage time x 



2. Limit Td 







In this limit, the first exponential term between the 
large curly brackets in Eq. (|40p. which contains the ex- 
pression (s + 1/Td), vanishes. The reason is that the 
real part of the exponent tends to —oo as Td tends to 
0. Therefore, the linear correction term of the density 
function, in the Laplace domain, simplifies to 



4f(.) 



2D 



r (a;th 
exp I 



2D 



}.(53) 



As in the previous limit, the linear term in the ex- 
pansion of the moments can be obtained directly from 
Eq. (|53|) and coincide with the corresponding case in [9| . 

The simplified expression obtained, Eq. (f53| . is ana- 
lytically tractable and the inverse Laplace transform is 
easy to compute [l^ . In the temporal domain, the linear 
correction is 



► 0) 



[Xthr 



xq) 



VAttDt^ 

1 



exp 



2Dt 

Mr) 
2Dt 



{ 



[(a;thr 



xa) 



^^rf 



{Xthi 

{Xthr - Xq) - flT - 



ADt 
2Dt 

{xthr - Xo) 

2Dt 



{xthr - Xo) 



.(54) 



In Fig. [5l^b) we can observe this limit expression as a 
function of t. Additionally, simulations based on differ- 
ent Td values show that the limit is reached for sufficiently 
small values. 



FIG. 5: (Color online) Limit behaviors for the linear correc- 
tion function (^i(r). (a) In the limit Td — >■ oo the product 
Td ■ <^i(t) does not depend on Td [solid (black) line], Eq. (|50[1 . 
The empirical product is shown for different Td values (col- 
ored stair- like lines) . As Td increases, the empirical histogram 
(obtained as in the insets in Fig. [T] multiplied by the corre- 
sponding Td) coincides with the analytical expression, (b) In 
the limit Td — the density function (j>i{j) does not depend 
on Td [solid (black) line], Eq. (|54p . The empirical function 
is shown for different Td values (colored stair- like lines). As 
Td decreases, the empirical histogram tends to the analytical 
expression. Parameters of the simulation as in Fig. [Tl with 
e = —0.1 and different Td values. 



From Eq.l54| we can compute the FPT density function 
for this limit which, up to order 1, reads 



[r) 



2Dt 



{xthr - Xo) - AiT - 



7 ) Mr)- 

■thr - Xq)\ J 

(55) 

The limit Td — >■ also enables a physical interpreta- 
tion. Since exp[— (i — to)/Td]/Td — 5(io) as Td — >■ 0, from 
Eq. ([1]) it is easy to see that after a differential time from 
the initial time Iq + dt, the system has moved to the po- 
sition xq H-e, and thereafter, the dynamics is governed by 
a constant drift /i. In this case, the FPT density function 
is (I)o{t), Eq. (|151) . with the initial position modified: 



[{xthi -Xp-e)- i.itY 
ADt 



(t) - (^thr - Xq - e) 

VAttDt^ 

(56) 

The expansion of Eq. (|56p around e — up to order 1 
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also coincides with Eq. ([55]) . 



IV. DISCUSSION AND CONCLUDING 
REMARKS 

In the present work we have analyzed the survival 
probability and the FPT problem of a Wiener process 
driven by an exponential time-dependent term superim- 
posed to a constant drift, Eq. ([T]), in the presence of an 
absorbing fixed boundary. We first focus on the survival 
probability in the region of interest and derive the time- 
inhomogeneous diffusion equation governing it, Eq. 
in the framework of the backward FP formalism. We 
propose a solution given by an expansion in terms of the 
intensity of the exponential drift, Eq. and derive the 
associated equations and (boundary and initial) condi- 
tions to solve each term, Eqs. dS]) to ([TT|) . Interestingly, 
the resulting equations are recurrent and easy to solve 
via a Laplace transformation. We explicitly solve up to 
the second-order term in the expansion, in the Laplace 
domain, and give some remarks about higher-order terms 
(see corresponding subsections). In particular, we show 
that each term exists, and therefore, the expansion we 
proposed is justified and constitutes the exact solution. 
Moreover, when the solution is set to the initial condi- 
tions of the problem, the probability depends exclusively 
on the time elapsed from the initial time, as expected 
from physical considerations, Eq. (I^Hl) . The FPT density 
function is obtained in terms of the survival probabil- 
ity, and we show that the expansion is preserved in this 
function and its moments, Eqs. ([M]) and (|42l) . Since the 
solution of each term is easily obtained in Laplace domain 
and the inverse transform is not always available, we re- 
view some related properties that can be calculated from 
them: the moments of the density function and the spec- 
tral density of an associated renewal process or "spike 
train" . 

In the second part of this work, we focus on the com- 
parison between the explicit results we have obtained and 
numerical simulations. Since truncation of the series re- 
sults in an approximate solution, we mainly focus on the 
first-order expansion. This linear regime coincides with 
a perturbation scenario. As shown in Fig. [Ija), the first- 
order term in the expansion of the FPT density com- 
pletely defines a slightly perturbed system, whereas a 
higher intensity of the exponential time-dependent drift 
facilitates higher-order effects [Fig. [T](b)] . Linear expan- 
sion of the first two moments of the density function 
reproduces numerical results accurately, except for low 
time constants of the exponential drift (Figs. [H and [3]). 
For a small exponential drift, we calculate the spectral 
density of the resulting spike train and observe that ana- 
lytical results coincide extremely well [Fig. |4l[b)]. More- 
over, the change in the spectral properties due to a small 
perturbation can be mostly ascribed to the change in 



the mean FPT [Fig. HKc)]. Finally, we derive the be- 
havior of the linear expansion of the FPT density func- 
tion, in the limit of negligible as well as extremely large 
time scales, for the time constant of the exponential drift. 
These limit expressions are inverse Laplace transformed 
and compared with exact results obtained from physical 
considerations, valid strictly in the corresponding limit 
(Fig.©. 

The process considered in this work naturally arises in 
neuroscience, but it is not restricted to this field (e.g., 
consider the motion of a charged particle in an exponen- 
tially decaying electrical field). In the context of neuro- 
science, the behavior of stochastic spiking neurons with 
an adaptation current can be described by stochastic pro- 
cesses with an exponentially decaying temporal term [lq|. 
For spiking neurons, the state variable x corresponds to 
the membrane potential and its evolution is given by a 
Langevin equation, where the integration of an input cur- 
rent is performed until a threshold is reached. At this 
moment a spike is generated and the time corresponds 
to the FPT in the statistical description [H HI]. The 
Wiener process is the stochastic representation of the 
basic model in theoretical studies, namely, the perfect 
integrate-and-fire neuron. In addition to external signals, 
an adapting neuron integrates, between spikes, an expo- 
nential time-dependent current corresponding to specific 
ionic channels [26;]. Therefore, the statistical description 
of the model we have considered provides important mea- 
sures for analyzing further effects of adaptation on spik- 
ing neurons. 

Finally, we note some features about the methodology 
considered here. First, we derive the equation governing 
the transition probability for a specific temporal drift, 
Eq. ([SJ, and propose a solution in terms of a certain ex- 
pansion, Eq. ([5]). Both procedures can be used for any 
time-dependent drift. However, in general, the resulting 
equations would have a source term difficult to tackle 
analytically. Second, our starting model is the Wiener 
process with an exponentially decaying temporal drift. 
It is easy to check that an equivalent formulation can 
be performed to other one-dimensional processes (e.g., 
Ornstein-Uhlenbeck process) with the same kind of tem- 
poral drifts. In these cases, the only change with respect 
to the Wiener process is that the homogeneous parts of 
the differential equations we have obtained, on the left- 
hand side of Eqs. ((H]) and ©, are different. Obviously, 
the difficulty in solving these other cases depends on the 
system at hand. 
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